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Understanding magnetism and electron correlation in many unconventional superconductors is 
essential to explore mechanism of superconductivity. In this work, we perform a systematic numerical 
study of the magnetic and pair binding properties in recently discovered polycyclic aromatic 
hydrocarbon (PAH) superconductors including alkali-metal-doped picene, coronene, phenanthrene, 
and dibenzopentacene. The ^-electrons on the carbon atoms of a single molecule are modelled by the 
one-orbital Hubbard model, and the energy difference e between carbon atoms with and without 
hydrogen bonds is taking into account. We demonstrate that the spin polarized ground state is realized 
for charged molecules in the physical parameter regions, which provides a reasonable explanation of 
local spins observed in PAHs. In alkali-metal-doped dibenzopentacene, our results show that electron 
correlation may produce an effective attraction between electrons for the charged molecule with one or 
three added electrons. 



The recent discovery of superconductivity in the solid state of various alkali-metal- doped PAHs, including 
picene 1 , coronene 2 , phenanthrene 3 ' 4 , and 1,2:8, 9-dibenzopentacene 5 , has attracted a growing interest in the 
condensed matter physics community. One of the most intriguing properties is that the superconducting 
transition temperature T c increases dramatically with increasing the number of benzene rings, reaching up to 33.1 
iCas reported in alkali-metal-doped 1,2:8, 9-dibenzopentacene 5 , comparable to that of alkali-metal-doped C 60 6 " 8 . 
This behavior suggests that PAHs with longer chains of benzene rings may exhibit higher T c . Another interesting 
experimental finding is that the magnetic susceptibility in the normal state of superconducting compounds is 
considerable large, and shows the Curie-Weiss behavior as a function of temperature 1,3 , evidencing that there exist 
local spins in the charge doped systems. Currently, the physical origin of local spins is unclear. 

Theoretical studies mainly based on the first-principle calculations have been carried out to elucidate the 
superconducting mechanism in these PAH superconductors. There are three important findings: (l)The con- 
duction band group around the Fermi level comprises four bands, which are derived from the lowest unoccupied 
molecular orbital (LUMO) and LUMO + 1 9 " 11 . (2) The bandwidth W around the Fermi energy is smaller than the 
effective Coulomb interaction U e ff in alkali- doped picene, placing it in the strong correlation regime 1213 . (3) 
Intramolecular 14,15 , intermolecular, and intercalant 16 phonons are shown to have significant contribution to the 
electron-phonon interaction in alkali-doped picene. However, Kato et aV 7 found that the electron-phonon 
coupling constant decreases with an increase of the number of carbon atoms in phenanthrene edge-type hydro- 
carbons, leading to a decrease of T c with increasing size of phenanthrene edge-type hydrocarbons. This contra- 
dicts to the experimental results that T c increases from ~ 5K in phenanthrene to ~ 33K in 1, 2 : 8, 9- 
dibenzopentacene. Finding (3) indicates that the superconductivity in PAHs can not be explained by considering 
only the electron-phonon interaction. Finding (2) suggests that electronic correlations are required to be taken 
into account for a complete understanding of the physical properties of PAHs. 

To gain some insights on the effect of electronic correlations on the magnetic and pair binding properties of 
PAHs, we undertake a first step to study the electronic states on a single molecule, which provides an energy scale 
intermediate between the microscopic energetics at the level of carbon atom and the macroscopic scale of the 
molecular solids. The resulting many-body states form a basis for analyzing the low-energy physics at the length 
scale larger than the size of molecule 18 . The carbon atoms in PAHs are in a trigonal planar sp 2 hybridization 
configuration with a half- filled p n orbital perpendicular to the planar structure, which forms energy levels 
(including LUMO and LUMO + 1) around the Fermi energy. The minimal description of interacting /^-electrons 
on a single molecule is given by the one-orbital Hubbard model which incorporates an on-site Coulomb repulsion 
besides a simple Huckel picture for the one-electron part. The Hamiltonian reads, 
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The operators c\ a and c ia create and destroy a n- electron with spin a 
at the ith carbon atom of a single molecule, respectively. n io = c] G c io is 
the number operator for an electron with spin a at site i, and = 
»^+n,-j. The sum ((/) is over nearest-neighbor (NN) carbon atoms, 
and the primed sum runs over the carbon sites without connecting to 
hydrogen atoms. Here, U denotes the on-site Coulomb interaction 
and e the energy difference between carbon atoms with and without 
hydrogen bonds, which is due to different electron-negativity of 
carbon and hydrogen. Notice that the on-site Hubbard U is different 
from the effective Coulomb interaction U e ff 12 ' 13 , which stands for the 
energetic cost to doubly occupy a molecular orbital. Based on the 
estimated parameters 18 " 20 : t ~ 2.5 eV to 3.0 eV and U ~ 6 eV to 
10 eV, we mainly used U = 2t and 3t in our studies. Two accurate 
numerical methods, exact diagonalization (ED) 21 and constrained- 
path Monte Carlo (CPMC) 22 ' 23 , are employed to simulate the model 
Hamiltonian (1). In the folio wings, the results presented are obtained 
by the CPMC method, except where explicitly noted otherwise. 

Results 

Four hydrocarbon molecules investigated in our studies are dis- 
played in Fig. 1, which are defined as C 14 H 10 (a), C 30 H 18 (b), 
C22H14 — A (c), and C 22 H 14 — B (d). The former three correspond 
to the hydrocarbon molecules in phenanthrene, 1,2:8, 9-dibenzo- 
pentacene, and picene, respectively. C 2 2H 14 — B, which has a similar 
structure to C 30 H 18 with both linear and angular fusion of benzene 
rings, is degenerate with C 2 2H 14 — A from our first-principle calcu- 
lations, suggesting its possible existence in alkali-metal-doped 
picene. 

In the noninteracting limit, i.e. (7=0, the electronic spectrum can 
be easily obtained by diagonalizing the quadratic Hamiltonian in Eq. 
(1). The energy difference between LUMO and LUMO + 1, 
E(LUMO+l)-E(LUMO), is shown in Fig. 2 as a function of e for 
different molecules. It is obviously that the electronic spectrum 
depends sensitively on the molecule structure. For C 14 H 10 and 
Q2H14 — A with phenanthrene edge- type termination, the energy 
difference decreases monotonically with increasing e; while for 
C 30 H 18 and C 22 H 14 — B, it decreases first with increasing e, and then 
turns to increase with e beyond a certain e. The energy differences 
between LUMO and LUMO+1 in phenanthrene and picene from 
the first-principle calculations are estimated to be about 0.1 eV and 
0.07 eV (see Supplementary Information), respectively, which are 
consistent with the results for C 14 H 10 at e~0.8f and for C 22 H 14 — 
A at e~0.87f. This indicates that PAHs could be well described by 
the model Hamiltonian (1) with e>0.8t. 




C 2 2H I4 -A 



C::H|4-B 




Figure 2 | Energy difference between LUMO and LUMO+1 as a function 
of e for Ci 4 H 10 (solid line), C 2 2H 14 — A (dashed line), C 2 2H 14 — B (dotted 
line), and C 30 H 18 (dashed-dotted line). The gray arrow indicates the 
position where E(LUMO + 1) — E(LUMO) is approximately 0.1 eV for 
C 14 H 10 and C 2 2H 14 — A. 

Magnetic property. First, we show the spin phase diagram of 
different molecules in Fig. 3 at U = 2t. Here, we present the results 
for the neutral molecule (the total number ofp^-electrons N e equals 
the total number of carbon atoms N c ) and charged molecules with 1, 
2, 3, and 4 added electrons. The total spin S of the ground state is 
obtained by comparing the energies in different spin sectors. The 
molecules not shown in Fig. 3 always lie in the lowest spin states (S = 
0 for even N e and S = 1/2 for odd N e ). One can clearly see that for the 
two electrons added cases, the total spin S switches from 0 to 1 at a 
certain critical e for all molecules shown in Fig. 1. Moreover, for the 
molecule C 30 H 18 , the four electrons added case shows a transition 
from the high spin (S = 1) state to the low spin (S = 0) state as e is 
increased. At U = 3t, the phase diagram is similar to the one 
presented in Fig. 3. Besides, we find that C 30 H 18 with three 
electrons added changes from being in the low spin (S = 111) state 
to the high spin (S = 3/2) state at e = 0.97f. 

The spin polarized state for the charged molecules with two added 
electrons at large e provides a unified explanation of local spins 
observed for PAHs superconductors in the normal states, with an 
average number of added electrons per molecule n ave ~ 3 1 " 5 . In alkali- 
metal-doped PAHs, the added electrons come from alkali-metal 
atoms intercalated between the stacked molecules 9 " 11 . Therefore, 
local inhomogeneous intercalation of alkali-metal atoms in the 
molecular solid may induce a local nonuniform charge 
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Figure 1 | Sketch of aromatic hydrocarbon molecules with (a) 14 carbon 
atoms, (b) 30 carbon atoms, and (c) and (d) 22 carbon atoms. 



Figure 3 | Spin phase diagram for charged hydrocarbon molecules as 
displayed in Fig. 1 at U = It. The total spin S of the ground state is 
shown for a combination of molecule structure, charge, and e. The integer 
in the bracket denotes the number of electrons added to the neutral 
molecule. 
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Table I | ^dependence of the energy difference AE = E(S= 1) — 
E(S = 0) for the charged molecules with two added electrons at U 
= 2t and e = l.Of. Statistical errors are in the last digit and shown 
in the parentheses 



V/t 


AE(C 14 H 10 ) 


AE(C 22 H 14 -A) 


AE(C 22 H 14 - B) 


0.0 


-0.0452(6) 


-0.0074(8) 


-0.1073(7) 


0.1 


-0.0305(6) 


0.0009(8) 


-0.1041(9) 


0.2 


-0.0016(8) 


0.0121(9) 


-0.1017(9) 


0.3 


-0.0044(8) 


0.0257(9) 


-0.0965(9) 


0.4 


0.0134(9) 


0.043(1) 


-0.090(1) 


0.5 


0.0289(9) 


0.060(1) 


-0.086(1) 



distribution on different molecules. If two electrons happen to lie on 
the same molecule, the spin polarized state is favored, which con- 
tributes to the high magnetic susceptibility observed in experiments. 

Now the question arising is why the Curie-Weiss behavior of 
magnetic susceptibility is only observed when n ave is larger than 
two 1 ' 3 . This issue can be addressed by studying the effect of long 
range Coulomb interaction on the ground state of charged molecules 
with two added electrons. An inclusion of NN Coulomb interaction 
V niHj (V denotes the strength of screened NN Coulomb repul- 
sion) to the Hamiltonian (1) indicates that the energy difference 
between the high spin (S = 1) state and the low spin (S = 0) state 
increases with increasing V, as shown in Table I. One can observe 
that this energy difference becomes positive for V ^ OAt and V ^ 0. 1 1 
on C 14 H 10 and C 22 H 14 — A, respectively, suggesting that the low spin 
state is more stable when the NN Coulomb interaction is not all 
screened out. Considering the fact that electronic screening decreases 
with decreasing the metallicity of the molecular solid, it was expected 
that when n ave is not large enough, NN Coulomb interaction is not 
screened out, making the low spin state more stable than the spin 
polarized state. 

Pair binding property. In a strongly correlated electronic system, 
electron pairing may be induced either by a purely electronic 
mechanism or by the interplay of electron -electron interaction and 
electron-phonon coupling. The effect of electronic correlations on 
the superconductivity can be explored by examining the pair binding 
energy, which is defined as: A z = lE—E^— E i+1 with E t repre- 
senting the ground state energy for the charged molecule having i 
added electrons 18 . In this notation, a positive A,- implies that an effec- 
tive attraction between added electrons can be produced through a 
charge disproportionation process between the two molecules: (/) + 
(/) —> (i — 1) + (i + 1), meaning that it is energetically favorable to 
have two molecules with i — 1 and i + 1 added electrons each rather 
than two molecules with i added electrons each. We have calculated 
the pair binding energy A f (/ = 1 and 3) for all molecules as shown in 
Fig. 1 and the results show that it is always negative in the whole 
parameter regime for C 14 H 10 , C 22 H 14 — A and C 22 H 14 — B (see 
Supplementary Information). In Fig. 4, we show A x and A 3 for 
C 30 H 18 as a function of e. It is interesting to see that both A x and 
A 3 take positive values at U = 2t and U = 3t for e> l.Of, demon- 
strating that electron pair is favored in the C 30 H 18 molecular solid at 
n ave — i o and n ave = 3.0. The positive A x and A 3 actually represent 
the energy gain due to the formation of bounded S = 1 doubly- 
charged C 30 H 18 and bounded S = 0 quadruply- charged C 30 H 18 , 
respectively, when two electrons are added to the neutral or 
already doubly- charged molecules. A similar pairing process was 
suggested to be responsible for superconductivity seen in alkali- 
meta-doped C 60 , but it occurs for U > 3.3f 18 . 

The effective attraction between added electrons resulted from the 
positive pair binding energy may lead to various charge density- 
wave, magnetic, and superconducting states, depending on the 
details of interactions between electrons around the Fermi level 
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Figure 4 | Pair binding energies A x (a) and A 3 (b) versus e for C 30 -Ffi8. 
The on-site interaction U is indicated by the shape of symbol. 

(including the intermolecular interaction which has not been con- 
sidered in our present study) 24 . In the case that instability to super- 
conductivity is predominant, it is expected that the positive pair 
binding energy has a contribution to the formation of Cooper pairs. 
In the potassium-doped l,2:8,9-dibenzopentacene, superconduct- 
ivity was observed for 3.0 < n ave < 3.5 5 . Therefore, a positive A 3 at 
large epsilon suggests that electron correlation might have a contri- 
bution to the formation of Cooper pairs. The bounded S = 0 state for 
quadruply- charged C 30 H 18 is consistent with possible s-wave 
superconductivity in dibenzopentacene, of which the magnetic sus- 
ceptibility shows a similar temperature dependence to the one in 
barium-doped phenanthrene, where s-wave superconductivity has 
been demonstrated by specific-heat measurements 25 . At e = 1 .Of , our 
calculations with an inclusion of nearest-neighbor Coulomb inter- 
action V KiKj show that both A x and A 3 are increased by a finite 
V, indicating that the effective attractive interaction between elec- 
trons is robust against the long range Coulomb interaction. 

Discussion 

Our numerical calculations show that there exists strong magnetic 
instability in the charged molecules of PAHs, especially for the even 
number of added electrons, which provides a reasonable explanation 
of local spins observed in experiments. In alkali-metal- doped diben- 
zopentacene, electron correlation may produce an effective attrac- 
tion between added electrons through charge disproportionation in 
different C 30 H 18 molecules. This is not only useful for understanding 
the superconducting mechanism in alkali-metal-doped dibenzopen- 
tacene, but also encouraging for the search of higher T c PAH super- 
conductors with longer chains of benzene rings. Experimental test on 
the spin state and charge distribution in PAHs can be compared with 
our predictions and also is crucial for clarifying the magnetic and 
superconducting properties of charged PAHs. 

Methods 

To study the magnetic and pair binding properties of PAHs, we have mainly used the 
CPMC method 22 ' 23 to compute the ground state of interacting p n - electrons described 
in the model Hamiltonian (1). The basic strategy of the CPMC method is to project 
the ground- state wave function from an initial wave function (with correct quantum 
numbers corresponding to the ground state) by a branching random walk in an 
overcomplete space of constrained Slater determinants, which have positive overlaps 
with a known trial wave function. For more technique details we refer to Refs. [22, 23] . 
Since the total spin S of the ground state is unknown a priori, we first applied the 
CPMC method to project the lowest eigenstates of Hamiltonian (1) in different spin 
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sectors (S, S z ), and then defined the one with lowest energy as the ground state. In our 
calculations, the free-electron wave function, given by Eq. (1) at U = 0 and repre- 
sented by a single Slater determinant, was chosen as the initial and trial wave func- 
tions. This choice achieves the projection in a certain spin sector with 

S= \S Z \ = - |Nf — Ni |. Here, and N± denote the total numbers of electrons with 

spin up and spin down, respectively. 

To justify the accuracy of the CPMC program, in Fig. SI and Fig. S2 (see 
Supplementary Information) we compare the pair binding energy obtained from both 
ED and CPMC methods on Ci 4 H 10 at U = It and U = 3t. We observe that the CPMC 
data are in good agreement with the ED results and the systematic error induced by 
the constraint of Slater determinant is small. In Ref. [22], extensive benchmark 
calculations for the one-band Hubbard model defined on a square lattice also showed 
that in the weak and intermediate coupling regimes, the CPMC method can present 
us accurate energy and reliable correlations of the ground state. 
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